Phonon overlaps, polaron overlaps, and the effects of molecular size on quantum 

tunneling in polyacetylene 
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We provide a theory for the effects of polarons and phonons in mediating and suppressing the 
quantum tunneling of electrons into single molecules of conducting polymers, motivated by exper- 
iments on molecular quantum dots. The effects of both phonons and excitations of the polaron 
particle-in-a-box states are calculated. Using the Su-Schrieffer-Heeger (SSH) model of polyacety- 
lene, we calculate the suppression of ground-state to ground-state transitions and the position and 
strength of the phonon and polaron side-bands. 
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I. INTRODUCTION 

Molecular quantum dotsii^ii^i^ and conduction 
along one-dimensional molecules (particularly carbon 
nanotubes^'^) have both attracted much attention in the 
nanophysics experimental community. Electron tunnel- 
ing into molecular quantum dots is known to occur to- 
gether with the emission of particular phonon modes. 
The relative rates of emitting different phonons has 
been calculated^ using overlap integrals of the phonon 
wavefunction. Tunneling into one-dimensional polymeric 
molecules will have an additional complication: the elec- 
tron final state will be a polaron, dressed by a local 
screening cloud of phonons. Polaronic effects, although 
likely small for stiff carbon nanotubes, are important in 
other one-dimensional conducting polymers. 

Here, we calculate the suppression of tunneling due to 
polaron formation in the conducting polymer polyacety- 
lene, using the Su-Schrieffer-Heeger (SSH) model^. We 
calculate the phonon sidebands directly and provide the 
theory for calculating also the excitations into polaron- 
in-a-box excited states which should form a characteristic 
size-dependent series of excitation lines. 



II. POLYACETYLENE 

Polyacetylene^ is an organic, quasi one-dimensional 
semi-conducting polymer. Many experiments have been 
performed to characterize the physical properties of poly- 
acetylene such as its band gap, dimerization constant, 
lattice spacing, behavior under doping, and crystal struc- 
ture (packing). In its stable, dimerized form, trans- 
polyacetylene has a backbone consisting of alternating 
single and double bonds between adjacent CH groups. 
This structure and the degenerate ground states result- 
ing from changing the order of single and double bonding 
within the molecule accounts for many of the unusual 
properties of trans-polyacetylene such as the formation 
of solitons. Here, because we add only one electron, the 
ground state will form an ordinary polaron rather than 



solitons. Hence, the unusual properties of polyacetylene 
are not central to our analysis: the polaronic effects we 
examine will be shared by more mundane conducting 
polymers. 



A. The Hamiltonian 

Polyacetylene can be simply modeled by focusing on 
the behavior of the free carbon valence electron occupy- 
ing the pz orbital. The semi-conducting properties are 
thus due to polyacetylene's half-filled n band. 

Following Su, Schrieffer and Heeger-SiiS, the electronic 
band structure can be calculated using the tight-binding 
approximation for the electrons in the Pz orbitals. The 
unit cell will have a basis consisting of two CH groups 
with one pz electron per group. The simple model is 
valid where the 7r-band width is large compared to the 
band gap, and where the short-range electron-electron in- 
teractions are small compared to the valence band width. 

The SSH Hamiltonian for the lattice model of poly- 
acetylene is: 



TC — Tie 



n. 



el—ph 



n 



ph 



(1) 



where there are contributions from the electrons in the 
TT-band, the 7r-band electron interaction with phonons, 
and the phonons. 

The tight-binding Hamiltonian written in bra-ket no- 
tation gives the electronic contribution: 

Hei = - ^ toiln + 1, s) (n, s\ + \n, s){n+l, s\] (2) 



where to is the hopping matrix element connecting near- 
est neighbors, spin s = ±i, and n denotes the index of 
the CH group for which |n, ±i) are the pz orbital states. 
The hopping matrix element in the absence of dimeriza- 
tion is represented by io- 

The electron-phonon interaction is: 
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where a is the electron-phonon couphng constant and the 
electron-phonon couphng is expanded to first order about 
the undimcrizcd state assuming a smah dimerization. 
The lattice component is given by: 
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where the lattice coordinate m„ defines the horizontal 
displacement of the n*'' CH group away from its undis- 
placed lattice position at na with a the characteristic 
lattice length. M is the mass of the CH group and K is 
the spring constant for the cr-bonding backbone. For a 
dimerized chain, u„ = (— l)"uo, where uq is the dimeriza- 
tion constant. The constants to = 2.heV , K = 21ey/A, 
and a = 1.22A were chosen following Su, Schrieffer, and 
Heeger— . Writing down the Bloch wavefunctions for the 
band states and diagonalizing the electron Hamiltonian 
yields an expression for the energy in terms of the crys- 
tal wave vector k, the dimerization m, and the constants 
above. As the number of CH groups, N, tends to infinity, 
the sum in the expression for the electronic energy can 
be transformed into an integral; we solve it here for finite 
N. The SSH model spontaneously dimerizes to lower its 
energy, with uq « ±0.04A. 



III. FINITE LENGTH CHAINS 

Why are phonon sidebands seen for electronic transi- 
tions into molecular quantum dots when they have not 
been important in the larger mesoscale dots? We ad- 
dress this question by comparing the phonon overlaps 
for molecules of different sizes/lengths. To be specific, 
we consider finite chains of polyacetylene with a discrete 
number N of CH groups. For simplicity, we consider the 
case of periodic boundary conditions such that the first 
and last CH groups are linked. Given that we are look- 
ing at dimerized polyacetylene chains, N must be even 
to avoid imposing a soliton into the molecule. In ad- 
dition, the properties of small, circular SSH molecules 
are known to have a sharp alternation between lengths 
N — An and A^ = 4n -|- 2; we focus here on lengths which 
are multiples of four to avoid distraction by these finite- 
size oscillations. From the hessian, we can extract a 
dispersion relation for various lengths of chain. Figure Q 
shows the dispersion relations for charged polyacetylene 
chains of length N=28 and N=48. Note that the range 
that the frequencies fall in are about the same for both 
chains (0-32 tHz). However, their dispersion relation dif- 
fers. The acoustic branches (neighboring carbons moving 
in concert) are rather similar. The optical branch (alter- 
nating displacement directions for neighboring carbons) 
of the shorter chain exhibits a classic dispersion relation 
for a dimerized system^, but the longer chain shows a 
pronounced softening, foreshadowing the formation of a 
polaron for chains of A > 50. Notice in particular the 
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FIG. 1: Dispersion relation for charged polyacetylene chains 
of length N 



softening of the optical branch at the mode ko with one 
wavelength in the length of the chain; the polaron when 
it first forms near A^ = 52 has a roughly sinusoidal form 
with this wavelength (see Fig. 3). The zero-phonon emis- 
sion line for an electronic transition is suppressed by the 
square of the overlap integral Oo,o between the initial 
and final vibrational ground states. This is traditionally 
measured using the total (;- factor: 



G = -log(Poo) = -log[(Ooo)'] 



(5) 



where Poo — ^ is the probability that the molecule will 
remain in its ground vibrational state with the addition 
of an electron. 

The phonon sidebands and extra lines on the dl jdV 
curves arise from overlaps between the initial ground 
state and excited final states. The integrated strengths of 
these transitions (the sum of the heights for the vibration- 
induced steps on the I- V curve) is given by 1 — Pqo ■ Thus 
large G is a signature of the importance of phonon emis- 
sion. We first focus on the dimerization mode, changes of 
which dominate the lattice deformation for small CxHx 
rings in the SSH model. 



A. Dimerization mode only 

We start by restricting our attention purely to the 
change in dimerization uq in the SSH model as the N- 
monomer system changes from neutral to charged. To re- 
alize this simple picture, we treat the polyacetylene chain 
as a one-variable problem, optimizing only the dimer- 
ization constant uq- We find the minimum uq for neu- 
tral and charged chains of various lengths A^ by scan- 
ning the energy over different values of uq and constrain- 
ing the displacements of CH units to be of the form 

Un = (-l)"Uo. 
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FIG. 2: Chain length dependence of the phonon overlap G, 
considering the dimerization mode only. For large A'^, we find 
the overlap oc 1/A'^, as discussed in Appendix A. For A'^ > 50 
polarons form. 



Specifically, the dimerization of the initial neutral state 
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has an associated net Euclidean displacement dis- 



tance in A^-dimensional configuration space of ri = 
Ug WN. Similarly, r2 = Uq WN for the charged final 

state, so the net displacement is d = (uq — Mq )VN 
for our oscillator (A^ orthogonal coordinates all shifting 
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by Mg — Uq ) . The frequencies are given by the optical 
phonon frequency (where neighboring CH units move in 
opposing directions) in the neutral (wi) and charged {UJ2) 
states. The one-dimensional formula for the overlap is: 
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where the displacement d and the frequencies uji and W2 
are defined above. For example, this value of the overlap 
yields Oo,o = 0.468 and G = -log|Oo,oP = 1-52 for a 
twelve unit polyacetylene chain. 

Note that G falls off as -^ for large A^ in figure ^ . In 
appendix A, we consider broadly the finite-size effects of 
both optical phonons (considered in the above calcula- 
tion) and acoustic phonons. We will see quite generally 
that optical phonon overlaps for delocalized electrons go 
to one as A^ ^- 00 (G ~ 1/A^), while acoustic phonons in- 
duce overlaps which go to a constant for large A^. Hence, 
in the absence of polaronic distortions, optical phonon 
emission becomes less important for larger molecules, 
perhaps explaining why phonon sidebands first became 
significant when the dots became molecular. 



B. Beyond dimerization: polaron formation 

Especially in quasi one-dimensional systems, polaron 
formation upon injection of an electron becomes an im- 
portant consideration. For this section of the paper, we 
treat the polyacetylene as an N-dimensional problem, al- 
lowing us to study this polaron formation. 

Polarons are quasiparticles formed when an electron is 
self-trapped in a lattice deformation. Polarons are thus 
electrons dressed in phonons and affect the conductivity 
and other properties of a material. 

We find that polarons don't form for short chains (the 
polaron doesn't fit since its equilibrium length is less than 
than the chain length). They do form for longer chains 
and lead to a finite phonon overlap suppression and hence 
significant phonon sidebands during electron tunneling. 
Because of the presence of these polarons, these side- 
bands do not disappear as A^ -^- cxd. Similar important 
polaronic effects should be expected for other conduct- 
ing polymers and for 3D systems that form polarons with 
substantially suppressed electron mobility.-- 

(This is not to say that real, short polyacetylene rings 
do not distort when charged. Indeed, benzene (A^ — 6) 
undergoes a dynamic Jahn- Teller distortion when a hole 
is injected,"'^'^ into a shape which could be viewed as a 
soliton-antisoliton pair. The SSH model applied to ben- 
zene dimerizes, however, lifting the electronic degener- 
acy that drives the Jahn- Teller distortion. The quantum 
phonon fluctuations ignored in the SSH model (the res- 
onance between the two dimerized benzene states) are 
necessary to restore the hexagonal symmetry to neutral 
benzene, but are too small to destroy dimerization in 
long polyenes. Adding quantum fluctuations to the SSH 
model would treat charged benzene as a resonance be- 
tween two dimerized states; the more complete model 
(including bond angle rotations) treats charged benzene 
as a dynamical resonance between Jahn- Teller distorted 
states when phonon quantum fluctuations are added.) 

We first find the minimized geometry of the polyacety- 
lene chain with the polaron. The shape of the polaron 
is traditionally derived from a continuum theory, with 
the displacement field (j)n as a sum of two arctangents 
(reference [U, eqn 4.43): 



0n = (-l)"Un 
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where ^ is the width of the polaron, b — .623^, a is the 
lattice constant, and 0„ is the 'staggered' dimerization, 
a more continuous variable. Values for these quantities 
from the literature^'* are ^ = 7.017 A, and ug ~ 0.04A 
for the charged molecule. Figure [3] shows this traditional 
continuum ansatz for the polaron shape. 

However, we found that this traditional form did not 
describe our relaxed configuration, which showed an ad- 
ditional distortion (Fig. 3]) . The polaron in equilibrium 
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we reproduce this distortion (Fig. (|4])). For long chains 
(length^ 1000), we find the polaron width parameters are 
^ — 9.5A and b — 4.8A, the polaron amplitude increased 
by a factor A ~ 1.3, and the lattice contraction due to 
the polaron is i? = — 0.036A spread over a distance ^2 — 
11. 3 A. (This corresponds to a shrinkage of the average 
lattice constant a by a fraction X — B/^ over the length 
^ of the polaron.) We will use the approximate distortion 
of equation ([8]) below to estimate the effective mass m* 
of the polaron. 



FIG. 3: Traditional polaron shape (equation (I?))). The shape 
for molecules of length ^ < L is closely approximated by one 
wavelength of a sinusoid; in this limit the polaron is an optical 
phonon with wave vector k = 2-k / L. 
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FIG. 4: When allowing the system to relax fully in the pres- 
ence of an extra electron, systems of larger than 50 atoms 
exhibit the polaron. In addition to the pinching in the mid- 
dle, the polaron state also shows a distortion of the lattice 
at the two ends of the molecule which is due to the periodic 
boundary conditions and the fact that the chain is contract- 
ing. The diamonds show the partially relaxed configuration 
using our functional form (eqn (|8])); the circles show the fully 
relaxed configuration minimizing A'^ coordinates. 



has a slightly contracted lattice constant, which pulls 
in the surrounding atoms. By relaxing the total energy 
within the following functional form 



(-l)"«o 1 + 



(A/V2)( 



tanh 




— tanh 



C. Phonon overlaps and polarons 

The general formalism for calculating phonon over- 
laps and emission probabilities is discussed in reference!^. 
There are certain subtle modifications needed here in or- 
der to deal with the polaron. 

When the initial or final electronic state of the sys- 
tem forms a polaron, one must generalize our formulas to 
deal with the almost zero-energy translation mode of the 
polaron. The hessian for the charged state of the poly- 
acetylene chain for A^ > 50 has two eigenvalues which 
are near zero. One is common to all lengths and charge 
states and corresponds to the translational symmetry of 
the molecule; the associated eigenvector has all compo- 
nents equal, and this mode can just be factored out in 
the calculation (by singular value decomposition or by 
hand). The second mode is quite precisely given by the 
derivative of the polaronic deformation (equation ([7]), (O, 
or Fig. (HI) with respect to the center sq of the polaron, 
which we will call the polaron translation mode PT: 






d(j)n/dxo 



yJ2m {d(t),n/dxof 



(9) 



To linear order, this eigenmode translates the polaron 
sideways, without overall displacment of the underlying 
lattice of atoms. Translating the polaron sideways by 
one lattice constant brings the polaron to an equivalent 
configuration (with periodic boundary conditions), and 
hence does not change the energy. There is a small energy 
barrier between these equivalent configurations (some- 
times termed the Peierls barrier), but it is tiny: expo- 
nentially small as the polaron becomes large compared 
to the lattice spacin g^^d^ . To an excellent approxima- 
tion, the polaron is a free particle, confined to the length 
of the chain. 

Treating the polaron as a particle in a box, we need to 
find its ground-state wavefunction and excitation ener- 
gies. To a linear approximation, shifting the polaron by 



}PT 



adds 



6pT/\JY.ni {d(t),n/dxof j d(j)„ / dxQ to the dis- 



placement field, which changes the center of the polaron 



xo by SpT/yJ2rn (d'Pm/dxo) . The coordinate Spt as- 
sociated with the polaron translation eigenmode of equa- 
tion ([9]) has a kinetic energy mcH Sprp/2. Changing co- 
ordinates from SpT to a;o, the kinetic energy associated 
with the polaron translation is: 



^y^m {d4'm/dxo) = L^JmFJmcH, so the correspond- 
ing wave-function for the polaron translation mode is 



^n{5pT) = -rfX, 



. I 1TICH ik„SpTy/mcH/m* 



(13) 



where 



KpT = mcH Sprp/2 

= mcH ^ (dfprn/dxof I ±1/2 

= m*xl/2, (10) 



mcH^idcprn/dxQf (11) 



and the shift in the polaron center x — xq given the am- 
plitude 5pT of the translational eigenvector epT is: 



X- XQr^ 5pT\/mcH/m*. 



(12) 



The effective mass of the polaron is light. We calculate 
the effective mass for a chain of length N = 1000, sum- 
ming equation (|11|). (with discrete derivatives stepping 
by two) using the analytic expression (equation[5|) for the 
shape of the polaron. We find that the polaron mass is 
1.24 X l{)~^mcH ~ 377ie, somewhat larger than the value 
4.44 X \Q~^mcH ^ LOTwe found without allowing for 
the lattice contraction^^ (equation [7] and figure [3]) rather 
than the fully optimized distortion of figure ([4])). To en- 
sure that the Born-Oppenheimer approximation is still 
valid in our problem, we also calculate the band-mass of 
the electron^^ with the expression m* = Ti /^p-, evalu- 



determined to make the square of the wavefunction in- 
tegrated over the size of the box equal to one. The fi- 
nal ground-state wavefunction \1/q „ for the charged poly- 
acetylene molecule with a polaron is a product of Gaus- 
sians for all of the other normal modes of the molecule, 
times this one-dimensional wave function, where the first 
subscript of ^o.n will be used to denote the number and 
type of phonons emitted, and the second will denote the 
energy level of the polaronic translation mode. 

It is convenient (and necessary for molecules with more 
than one type of atom) to change coordinates from posi- 
tions X to y = M^/-^(x — ri), where ri is the equilibrium 
initial configuration. Treating the mass in the SSH model 
as a CH mass, and adopting the notation of reference^, 
the final-state wavefunction becomes: 
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ated at fc = f and taking the derivative of the hopping ^'^^ eigenvalue near zero), fc = y 



2a 

integral. The band mass for the electron is O.lSrng. This 
value still allows for the Born-Oppenheimer approxima- 
tion since the light electrons have sufficient time to relax 
as the polaron distortion fluctuates. The ground-state 
energy for the polaron center-of-mass coordinate is zero 
with periodic boundary conditions. For a chain of this 
length, the energy h k^ /2m* for the first excited state 
k = 2tt/L is 6.29 tHz and the second excited state has 
energy 56.61 tHz, both in the experimentally interesting 
range near the phonon frequencies of Fig. (fT|). 

The perhaps more realistic fixed boundary conditions 
with ground state fc = f- and excited state k = ^ yields 
1.57 THz and 6.29 THz respectively. For electrons tun- 
neling into the edges of systems with non-periodic bound- 
ary conditions, there could also be a notable dependence 
of the electronic matrix element on the polaronic final 
state. 

To calculate the transition probability into phonon and 
polaron excited states, we need the polaronic wave func- 
tions. The one-dimensional wave- function for the po- 
laron center of mass is just ^'(s^o) = -jf Gxp{iknXo) , 
since we are using periodic boundary conditions. The 
corresponding wave-function in the original phonon- 
mode basis i/>{5pt) is confined to a box of size 



where J7 = M~2KM2 is the square of the frequency 
matrix (and K is the spring constant matrix) , Ca are the 
normal modes of ^2 with eigenvalues Wq {spt being one 

A = M^(ri-r2) 



is the net atomic displacement due to the transition, the 



primed k' = M 
wavevector, and 
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ipT is the rescaled polaron 



(15) 




a=3 



is the normalization of the Gaussian wavefunction with 
the zero-frequency polaron translation mode removed. 

The overlap between the uncharged initial state and 
the charged polaron-containing final state is then 
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(16) 



where the first subscript on Oo;o,n represents the initial 
(ground) state, the second subscript denotes the phonons 
emitted in the final state, and the third is the excitation 
level of the polaron center-of-mass wavefunction in the 
final state. 



G for dimerization-only and polaron calculations 
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FIG. 5: G for the case with polaron and without. 



We rewrite^ this overlap integral as a single Gaussian 
with quadratic form VI = {Vti + ri2)/2, centered at A = 
\n-^Q.2 -A: 
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We begin by considering the ground-state to ground- 
state overlap Oo;o,0; leaving the polaron in its ground 
state k„ = 0. The integral in the last equation of p?| is 
, and 



just -^ 

Oo;0,0 
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The total g-factor is defined as the suppression of the 
step height in the I-V curve due to the ground-state to 
ground-state overlap of the atomic wavefunctions, identi- 
cal to equation ([5]) save for the extra polaron component 
in the final wavefunction: 



to the larger shift in the geometry, the overlap becomes 
small and the total g factor becomes large in comparison 
to the situation without a polaron. 

Figure ([5]) is a plot of the g factors in the dimerization- 
mode-only calculation where no polaron is allowed, com- 
pared to the N-dimensional calculation where a polaron 
is free to form. The graphs coincide for the most part 
for N < 50, where no polaron forms. (The simple calcu- 
lation involving only the dimerization mode is approxi- 
mate even in the absence of a polaron in the final state 
because the other phonon modes can emit even numbers 
of phonons due to their frequency shifts, even though the 
configurational shift A has no component along them.) 
The polaron which forms for iV > 50 is associated with 
a fixed phonon overlap integral of about 4.5. 

It is natural to ask whether the SSH-model calcula- 
tions in this paper can be fleshed out into a quantitative 
picture, using modern density functional theory (DFT) 
calculations for real polyacetylene molecules (as has been 
done, for example, for C72 and C140 bi-fuUerenesi^) . 
Early DFT calculations on polyacetylene found no dimer- 
ization after careful convergence in fc-points^^; more re- 
cent calculations have shown that the dimerization de- 
pends also upon the basis set and the particular den- 
sity functional variant used^S. Charged polyacetylene has 
been studied^ only for effectively rather short systems 
(one charge per 30 carbons), and while the total length 
change was studied no polaronic or soliton deformations 
were described. In our initial explorationsi^ we were able 
to explore molecules of length A^ = 20 using a basis set 
which was lacking the diffuse orbitals needed to correctly 
describe the charged state of this conjugated molecule. 
These calculations showed a distortion in the charged 
state primarily associated with the overall length of the 
molecule, with little effect on the dimerization. 



IV. CONCLUSION 

We have calculated the phonon overlaps for polyacety- 
lene using the tight-binding SSH model. We examined 
the suppression of the main transition (electronic) by 
both the phonons and the polarons, a quasiparticle that 
is important in quasi-one-dimensional systems. The ef- 
fect of the polaron is to increase the bare G for chains 
that would otherwise have a G that tends to zero as TV in- 
creases. Under certain conditions (when the chain length 
A^ > 50 and the polaron-in-a-box energy is comparable to 
the phonon energies), the presence of polarons modulate 
the bare phonon to phonon transition via the overlap in 
the same way that phonons modulate the bare electronic 
transitions. 



G=-log[|Oo;0,( 



(19) 



The presence of the polaron dramatically alters the over- 
lap between the ground state vibrational wavefunction 
of the neutral polyacetylene and the ground state vibra- 
tional wavefunction of the charged polyacetylene. Due 
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d^x 



p ■ ^{5/N) 



(22) 



where n labels the lattice sites. The only contributions 
are from longitudinal phonons; picking x, along k for a 
system of size L (with N = [L/a)^) we have: 



Appendix A 

We can get a qualitative understanding of why phonon 
emission is more important for smaller quantum dots 
by considering how adding an electron to a system 
with N atoms changes the bond lengths and hence the 
phonon overlap. There are two qualitative cases to 
consider — the short-wavelength optical phonons and the 
long-wavelength acoustic phonons. In the context of 
polyacetylene, adding a charge to a molecule with TV car- 
bon atoms (ignoring polaron formation) will change both 
the dimerization amplitude (an optical excitation, con- 
sidered in this manuscript) and the overall length of the 
molecular (an acoustic phonon). Charge-induced changes 
in the length of irans-polyacetylene have been studied us- 
ing density functional theories2£. 

For a single phonon mode (and no frequency shift dur- 
ing the transition), G — g is given by the square of the 
atomic displacement x between initial and final equilib- 
rium states divided by four times the square of the zero- 
point fluctuations x^: 



G = 



X 



4x0^ 



2h 



(20) 



When an extra electron is added to a quantum dot with N 
electrons, it usually delocalizes over ~ N bonds. There- 
fore, its effect on the length of each bond will be diluted, 
of order (Sa/N) where a is the interatomic spacing and 
(5 is a constant of order one. 

The optical case is particularly straightforward: N in- 
dependent atomic coordinates each shift by Sa/N, so the 



lidean distance x'^ = N{Sa/N)'^ = 


a^S^N. 


Hence 


g-factor becomes: 






x^ 52 1 

4x2 ^ ^• 




(21) 



Thus implying that optical phonon emission becomes 
smaller for larger quantum dots. We see this also in the 
explicit SSH calculation summarized in figure (P|). 

For acoustic phonons in a D-dimensional quantum 
dot, we must decompose the displacement field w(x) = 
{6/N)x of a dot uniformly rescaled in size by a factor 6/N 
into the orthonormal phonon modes ik — pe^^^"/^/N 
where p is the polarization vector (transverse or longi- 
tudinal). Each mode is excited by an amount qk where 
u = J2 Qk^k with 



Qk = efe • u 



E- 



■p-^n{5/N) 



r.L/2 



qk 



6 1 



dxixie""'-"^ L 



-D-l 



i-L/2 



N\k\ 



(23) 



The total g-factor is the sum of the factors gk for each 
normal mode (ignoring for simplicity the frequency shifts 
due to the extra electron), where we approximate Wk ~ 
ck with c the speed of sound: 



G 



E 



mujk 2 



E 



mck S^ 1 
^^7Vfc2' 



(24) 



We can approximate this sum as an integral, where the 
density of points in k-space is (2:^)^ , and make a spher- 
ical approximation to the Brillouin zone and the lower 



cutoff j^ < 



k\ < 



G 



zvr 



d^fc- 



N \ 
For D > 2, we have 



?71C\ 



2hk 

l/a 
1/L 



.D-2 



dk. 



(25) 



G 



L 



^52 .^. 



L 



N \2hJ \a 

D 



D-l 



N 2% 



= 5' 



2n 



and for D = 1 we find: 



G = 5^^a{\og{L/a)) 



(26) 



(27) 



We will find in 13 = 1 that the electron forms a polaron 
and is hence not delocalized, so this logarithmic diver- 
gence is bypassed in practice by polaron formation. 



The constant 5"^^ 



[aS) 



2 rauj£t 
2% 



{asy 



where the 



Debye frequency cod ~ c/a. Hence the acoustic phonons 
suppress the quantum tunneling by a finite amount . i , 
even for large quantum dots. 

Thus the emission of optical phonons during electronic 
transitions becomes unimportant for large systems, while 
the emission of acoustic phonons will continue to suppress 
the transition by a finite overlap integral even as A^ — > 
00. 
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Appendix B 

We now turn to transitions to final states with polarons 
in excited translational modes (with polaron wavevectors 
kn), and with zero and one phonon emitted. 

The integral in the last equation of (fT7|) is just the 
Fourier component of a multidimensional Gaussian at 
wavevector kj^ — ,^T", epT- Just as the Fourier trans- 



form of 



■ exp 



{-&) - exp(-^ 



so m many 



dimensions the Fourier transform of exp{~x'^ Ax) 
{l/y/det{A/TT)) exp(-k'^^-ik/4). Hence 



Oo;0,n 



dot 



cxp ( -j^K'^n-'K ) e'^""^ 



„ ^1^2 ^-lA^nA^-j,A^Q.A 






,e<^exp 



4^k„ i2 k„ 



Oo;o,o e^""'-^ exp L^k^^Q-^k; ) (28) 



This k-dependence is due to a kind of spatial confine- 
ment of the polaron center of mass in the overlap integral 
calculation. Indeed, when we changed coordinates from 
the polaron center of mass x to the translation eigen- 
mode SpTi we made an approximation that the polaron 
did not move substantially from its reference position 
xq. Polaron translations are given by adding multiples 
of dcfim/dxQ only for shifts x — ccq small compared to the 
polaron width ^. (A shift by many ^ would look like the 
sum of a separated polaron and an antipolaron, not the 
derivative of the polaron.) This approximation is valid 
for our overlap calculations so long as the 'zero-point' 
Gaussian fluctuations ^o of our overlap integral along x 
are small compared to ^. The RMS fluctuations in the 
polaron center-of-mass coordinate in the overlap integral 
are given roughly by 



{{x-xof)^ 



1 



- —A^lT) = —z{{y'^pTf 



m 



ra 



^^ j d^y{y ■ epT) 



1 T. 



xexp(--y'^Qy). 



(29) 



where the expectation is taken in a harmonic potential 
whose frequency matrix is fi, and we ignore the displace- 
ment A (which does not shift the center of the polaron) . 
But one can verify that J dy{y ■ v)^ exp(— y'^Cy) — 
iv^C~-^v (say, by going to normal modes of C, as in 
deriving equation [551 below), so the zero-point polaron 
fluctuations in our overlap integral, ^q, satisfy: 



eo 



{ix-xof)n 



2m* 



eprp 



n 



ipT- 



(30) 



For our polaron with N = 64, ^o = 16.56A« ^, so our 
linearized approximation is only approximate. (Compar- 
ing, for example, u{xq + ^o) ^ ""(^^o) in the SSH model 
with Mo -l- Co 1) shows a similar distortion, but with the 
former centered at xo + Co/2). 

Returning to our overlap integral in equation (|28p , and 
noting again that kj^ = ,^T", £pt, we see that the prob- 
ability of ending in a polaron excited state n with no 
phonons is writeable in terms of this zero-point polaron 
center fluctuation: 



Po:0,n — |Oo:0,n| — |Oo;0,n| 6 

— I '^0:0,0 1 e "■ L ' . 






Ip 



(31) 



And so the g factor associated with this transition is 
then: 



5o;0,„ = e--^"^(^)^=e-'=^ 



2 ^2 




(32) 



where fc„ — ^^. The initial step in the IV curve, where 
the polyacetylene molecule is left in its (charged) ground 
state is thus followed by a series of smaller steps at ener- 
gies -n — ? and relative heights e~^^^° . 

Finally, we turn to the case of one phonon emis- 
sion combined with excitation of the polaron translation 
mode: 



Oo;i.,„ = / rf^y*^'^*(y)*S,„(y) 



N^ 



=»<-y 






d^ye'k;.y~.^^(2)g-iy-Qy 



} 



la^ha. 



1^22 , 



V2WaA((A 




(33) 



The integral in equation (|33|) can be performed by 
decomposing y into the normal modes pp of il with 



{} 



E 



P/5-eL^^ / d2//3 2//3e-^""'-'^e 



dyp.e- 



Vli 






(34) 



Now j dxe^^^e 

ik 



IKCC ^-y — jt.3? 



re iA 



~ V A 

^ fc2 

,^ Y ^e 4A . And the overlap becomes 



and J dxe^'^^xe 



{} 



?(2) 






iK 



P(i 



2{uo/h) 

1 



n 

/3' 



ujp> /n 






(35) 



So, 

Oo;la,n 



,ih. 



giving us a g-factor for the combined excitation of a 
phonon and a polaron of 



Jk' A 



yZW^(^k;^fi^i.t) - (A - A) . e^)e 

Xe"4k„ !^^ k„ 1- 2 p-n^ S:2Ag-2jj-A O2A 



5la,' 



,i?i. 



N^\/% 



V^^JJhq^In-^e^ - (A - A) . 4^)) 
V2:W^(|kf i^-^ef ) - (A - A) . .t))Oo;o 



2m 



2 t2 l'2u)a 



-KC ( r^ (( 



,-(2; 



(37) 



huj 



2m 



-k ?^ 0-lf(2) 






2wo 



(A - A) . g 



?(2 



The first term is the zero-polaron phonon sideband 
suppressed by the polaronic overlap in the same way that 
in equation P^ the zero-phonon sideband is supressed 
by the phonon overlap. The second term is due to the 
change in the shape of the wavefunction from the pola- 
(36)ronic excitation. 
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